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—4- , We attempt to put constraints on different cosmological and biasing 

[--««. ■ models by combining the recent clustering results of X-ray sources in 

the local (z < 0.1) and distant universe (z ~ 1). To this end we compare 
the measured angular correlation function for bright [1] and faint [20] 
RC-SATX-ray sources respectively with those expected in three spatially 
y-H ' flat cosmological models. Taking into account the different functional 

q. forms of the bias evolution, we find that there are two cosmological 

models which match well the data. In particular, low-f2 cosmological 

models (Qa = 0.7) that contain either (i) high <r™ ass = 1.13 or (ii) low 

^_ > ■ <jg nass = 0.9 respectively with different bias behaviour, best reproduce 

•*» ' the AGN clustering results. While rCDM models with different bias 

behaviour are ruled out at a high significance level. 
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Introduction 

The study of the distribution of matter on large scales, based on dif- 
ferent extragalactic objects, provides important constraints on models of 
cosmic structure formation. In particular Active Galactic Nuclei (AGN) 
can be detected up to very high redshifts and therefore provide infor- 
mation on how the X-ray selected sources trace the underlying mass 
distribution as well as the evolution of large scale structure. 

However, a serious problem here is how the luminous matter traces 
the underlying mass distribution. Many authors have claimed that the 
large scale clustering pattern of different mass tracers (galaxies or clus- 
ters) is characterized by a bias picture [8]. In this framework, biasing is 



assumed to be statistical in nature; galaxies and clusters are identified 
as high peaks of an underlying, initially Gaussian, random density field. 
Biasing of galaxies with respect to the dark matter distribution was also 
found to be an essential ingredient of CDM models of galaxy formation 
in order to reproduce the observed galaxy distribution. Furthermore, 
different studies have shown that the bias factor, b(z), is a monotoni- 
cally increasing function of redshift [19]. For example, Steidel et al. [17] 
confirmed that the Lyman-break galaxies are very strongly biased trac- 
ers of mass and they found that b(z = 3.4) ~ 6,4, 2, for SCDM, ACDM 
(Q = 0.3) and OCDM (fi = 0.2), respectively. 

Studies based on the traditional indicators of clustering, like the two 
point correlation function, have been utilized in order to describe the 
AGN clustering properties. Our knowledge regarding the AGN clus- 
tering comes mostly from optical surveys for QSO's ([16] and reference 
therein). It has been established, that QSO's have a clustering length 
of r ~ 5.4 ± l.l/j -1 Mpc. 

Vikhlinin & Forman [20] studied the angular clustering properties us- 
ing a set of deep ROSAT observations. Carrera et al. [5] combined 
two soft X-ray surveys (235 AGN), found a spatial correlation length 
1.5/t -1 Mpc < r < 5.5/j -1 Mpc, depending on the adopted model of clus- 
tering evolution. Recently, Akylas et al. [1] using 2096 sources detected 
in the ROSAT All Sky Survey Bright Source Catalogue (RASSBSC), de- 
rived the AGN angular correlation function in the nearby Universe and 
utilizing Limber's equation obtained r Q = 6.7 ± 1.0/j -1 Mpc, assuming 
comoving clustering evolution. 

Here, we present the standard theoretical approach to estimate the 
angular correlation function w(9), using different models for the bias 
evolution in different spatially flat cosmological models. Comparing the 
latter with observational results, we attempt to put constraints on the 
different cosmological and bias models 

1. The Integral Equation of Clustering 

For the purpose of this study we utilize the relation between the an- 
gular w(0) and spatial £(r, z) two point correlation functions ([10] and 
references theirin). As it is well known, this connection can be done 
using the Limber equation [14]. For example, in the case of a spatially 
flat Universe (Q + Q\ = 1), the Limber equation can be written as 



,„ m _ o^Q Jo" N*(z)E(z)dz Jj? g(r, z)du 
{) ~ ~ US°N(z)dz]* • (1) 
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where x is the comoving coordinate related to the redshift through 

x = jrJ Et) ' (2) 

with 

E(z) = [0,(1 + zf + 0\] 1/2 (3) 

[14]. The mean surface density, A/", on a survey of solid angle Q s is: 

A/" = / x 2 4>{x)dx = — / iV(z)dz , (4) 

JO S2 S JO 

where iV (z) is the number of objects in the given survey within the shell 
(z, z + dz) and (j)(x) is the selection function (the probability that a 
source at a distance x is detected in the survey): 

/CO 
*(L x ,z)dL . (5) 

In this work we used a luminosity function $>(L X , z) of the form assumed 
by Boyle et al. [4]. 

The physical separation between two sources, separated by an angle 
9 considering a small angle approximation is given by: 

Extending this picture, we quantify the evolution of clustering with 
epoch presenting the spatial correlation function of the X-ray sources 
as 

Z(r,z)=S maBB (r)R(z) , (7) 

with R(z) = D 2 (z)b 2 (z), where D(z) is the linear growth rate of clus- 
tering (cf. Peebles 1993) being given by 

D{z) = ^^L -ww dy (8) 

and finally b(z) is the evolution of bias. 

2. Bias Evolution 

The concept of biasing between different classes of extragalactic ob- 
jects and the background matter distribution was put forward by Kaiser 
[8] in order to explain the higher amplitude of the 2-point correlation 
function of clusters of galaxies with respect to that of galaxies them- 
selves. Therefore, we shortly describe here some of the bias evolution 
models in order to introduce them to our analysis. 



2.1. No evolution of bias (BO) 

This model considers constant bias at all epochs: 

/ r \ 7 / 2 
*(*) = 6 ° - &?* [-^ ) (9) 

where r° p = 5.4/i _1 Mpc is the correlation length in comoving coordi- 
nates estimated by the APM correlation function [9] and r Q is the cor- 
responding length for the X-ray sources. Finally, 6° p = l/cjg iass is the 
present bias of optical galaxies relative to the distribution of mass and 
^.mass j g ^e mass rms fluctuations in sphere of radius 8h _1 Mpc. Using 
the above ideas, if ones assumes that r > r° p , then it is quite obvious 
that the X-ray sources are indeed more biased with respect to optical 
galaxies by the factor (r /ro P ) 7 ' 2 . In case that b(z) = 1, we have the so 
called non-bias model (from now on we consider b(z) = 1 as a (BO) bias 
model). 

2.2. Test Particle or Galaxy Conserving Bias 

(Bl) 

This model, proposed by many authors ([19] and references theirin), 
predicts the evolution of bias, independent of the mass and the origin 
of halos, assuming only that the test particles fluctuation field is related 
proportionally to that of the underlying mass. Thus, the bias factor as 
a function of redshift can be written: 

>w- 1+ %r ■ (10> 

where b Q is the bias factor at the present time. 

2.3. Merging Bias Model (B2) 

Mo & White [13] have developed a model for the evolution of the 
correlation bias, using the Press-Schechter formalism. Utilizing a similar 
formalism, Matarrese et al. [12] extended the Mo & White results to 
include the effects of different mass scales. In this case we have that 



( fee - 0. 41) 
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b(z) = 0.41 + ^JTs • ( n ) 



3. Cold Dark Matter (CDM) Cosmologies 

In this section, we present the cosmological models that we use in this 
work. For the power spectrum of our CDM models, we consider P(k) ~ 
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k n T 2 (k) with scale-invariant (n = 1) primeval inflationary fluctuations. 
We utilize the transfer function parameterization as in Bardeen et al. [2], 
with the corrections given approximately by Sugiyama's [18] formula: 

„, ln(l + 2.34g) r , , 2 

T{k) = 2 34g [1 + 3 ' 89q + {WAq) + 

{5A6qf + (6.71g) 4 ]- 1/4 . 

with q = k/hT, where k = 2tt/\ is the wavenumber in units of h Mpc -1 
and r is the CDM shape parameter in units of (/i — 1 Mpc) -1 . In this 
analysis, we have taken into account three different cold dark matter 
models (CDM) in order to isolate the effects of different parameters on 
the X-ray sources clustering predictions. 

The rCDM (fi = 1, h = 0.5, a^ 3 = 0.55, b = 2.18) and ACDM 
(0\ = 0.7, h = 0.65, erf™ 88 = 0.9, b = 1.33) [11] are approximately 
COBE normalized and the latter cosmological model is consistent with 
the results from Type la supernova [15]. In the same framework, the 
rCDM and ACDM models have T ~ 0.2, in approximate agreement 
with the shape parameter estimated from galaxy surveys [9] and they 
have fluctuation amplitude in 8 /i _1 Mpc scale, a™ 2188 , consistent with the 
cluster abundance, of 12188 = 0.55S1" - 6 [7]. 

Furthermore, in order to investigate cosmological models with high 
value of Og na88 , we include a new model named ACDM2 (Q\ = 0.7, 
h = 0.6, o-f™ 88 = 1.13, b = 1.06) [6]. The CTf ia88 value for the latter 
cosmological model is in good agreement with both cluster and the 4- 
years COBE data with a shape parameter T = 0.25. Therefore, it is 
quite obvious that two of our models (rCDM and ACDM) have the same 
power spectrum and geometry but different values of £l and cr™ 2188 , while 
the two spatially flat, low-density CDM models (ACDM and ACDM2) 
have different as and T respectively. It turns out that in ACDM and 
rCDM models the distribution of X-ray sources is "biased" relative to 
the distribution of mass; while the ACDM2 model is almost non "biased" 
while b at the present time is described by equation (9). 

In order to understand better the effects of AGN clustering, we present 
in figure 1 the quantity R(z) = D 2 {z)b 2 {z) as a function of redshift 
for the three cosmological models, utilizing at the same time different 
bias evolution. It is quite obvious that the behaviour of the function 
R(z) characterizes the clustering evolution with epoch. Figure 1, for 
example, clearly shows that the bias at high redshifts has different values 
in different cosmological models. In particular for the high o^ 88 low-Oo 
flat model (ACDM2) the distribution of X-ray sources is only weakly 
biased, as opposed to the strongly biased distribution in the rCDM 
cosmological model. 
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Figure 1. The function R(z) as a function of z, for different bias evolution models. 



Indeed the different functional forms oib(z), provide clustering models 
where: (i) AGN clustering is a decreasing function with redshift for 
(BO), (ii) AGN clustering is roughly constant for (Bl). However, the 
ACDM2-B1 model gives lower R(z) simply because the higher crg 1 ^ 
normalization largely removes the clustering difference between the two 
other flat cosmological models with low crg 1 ^ normalizations. In other 
words, the present bias value of the above model is almost ~ 1, which 
gives clustering behaviour similar to the ACDM2-B0. Finally, AGN 
clustering is a monotonically increasing function of redshift for (B2) and 
(B3) respectively ([3] for details). 

4. Application to the Data 

In Figure 2 (left) we compare the angular correlation function of a 
sample of 2096 sources with a total sky coverage of 4.9sr detected in the 
ROSAT All-Sky Survey Bright Source Catalogue [1] with that predicted 
in various flat cosmological models. Considering a two point angular 
correlation function of the form w(8) = {0/6 o ) l ~^ ', the above authors 
found O = 0.062°, 7 = 1.8 and spatial correlation length of r rj 6.5 ± 
1.0/i _1 Mpc and r ~ 6.7±1.0/t Mpc for stable and comoving clustering 
evolution respectively, similar to the optically selected AGN. Due to 
the fact that the above estimates have been focused on the local X- 
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Figure 2. Left: Comparison of the predicted angular correlation function for var- 
ious cosmological models with that of the local AGN distribution, estimated in [1]. 
Errorbars are determined by assuming Poisson statistics. The continuous dot-dash 
line represent the best fit to w x (9) derived by the above authors. Right: Comparison 
of the predicted angular correlation function for various models with that estimated 
in [20]. The continuous dot-dash line represent the best fit to w x (8) derived by the 
above authors. 



ray Universe, in this work we use complementary observational results 
from [20]. They have analysed a set of deep ROSAT observations with 
a total sky coverage of 40 deg 2 , in order to investigate the clustering 
properties of faint X-ray sources. Therefore, they claimed that the two 
point angular correlation function is well described by a power low with 
7 = 1.8. Thus, in Figure 2 (right) we plot their results and the estimated 
angular correlation function for all nine models. 

In order to quantify the differences between models and data, we per- 
form a standard x 2 test, for the bright (RASSBSC) and faint sources 
respectively, and we present the V^ 2 (see table 2), where we have been 
use N = B,F for bright (RASSBSC) and faint X-ray sources respec- 
tively. Comparing the statistical results for both (a) high (faint sources) 
and (b) low (bright) redshift regimes, we can point out that for the 
bright X-ray sources the only models that are excluded by the data, 
at a relatively high significance level, are ACDM-B2, rCDM-Bl and 
tCDM-B2. Interestingly, for the faint objects the excluded models are 
ACDM-B1, ACDM-B2, ACDM2-B0, ACDM2-B1, tCDM-B0, rCDM-Bl 
and rCDM-B2. The above differences between the two kind of popula- 



Table 1. \ 2 probabilities ("P^) of consistency between RASSBSC, faint ROSAT 
data and models. 

Comparison Pair 
RASSBSC - ACDM-BO 
RASSBSC - ACDM-B1 
RASSBSC - ACDM-B2 

RASSBSC - ACDM2-B0 
RASSBSC - ACDM2-B1 
RASSBSC - ACDM2-B2 



>x 

V B 2 

' >x 2 

0.81 

0.067 

0.019 

0.19 
0.42 
0.57 



Comparison Pair 
faint ROSAT- ACDM-BO 
faint ROSAT - ACDM-B1 
faint ROSAT - ACDM-B2 



V 2 

' >x 2 
0.037 

3.67 xlO" 4 

1.72 xlO" 6 



faint ROSAT - 
faint ROSAT - 
faint ROSAT - 



ACDM2-B0 
ACDM2-B1 
ACDM2-B2 



3.44 xlO 

3.21 xlO 

0.036 



RASSBSC - rCDM-BO 0.19 

RASSBSC - rCDM-Bl 2.88x 10" 
RASSBSC - tCDM-B2 3.00x10" 



faint ROSAT - 
faint ROSAT - 
faint ROSAT - 



rCDM-BO 
rCDM-Bl 
rCDM-B2 



4.64 xl0~ 
1.07x10" 
4.73 x lo- 



tions are to be expected simply because the cosmological evolution plays 
an important role on large scale structure clustering due to the fact that 
the high redshift objects are more biased tracers of the underlying mat- 
ter distribution with respect to the low redshift objects [17]. Also, from 
the faint X-ray sources results we would like to point out that there is 
not a single rCDM model that fits the data. 

If we make the reasonable assumption that there is no correlation 
between the two X-ray populations, mostly due to the large distances 
involved, we can consider the RASSBSC and faint ROSAT catalogues 
as being independent of each other. Under this assumption the previous 
statistical tests can be considered as independent (as indeed verified by 
a KS test). Therefore, the joint (overall) probability can be given by the 



following expression: P 



V B , 



1 >X 2 



1 - In (P. 



B 



V F - 
>X 2 >X 2 ' 



This overall 

statistical test proves that the ACDM-BO (P = 0.14) and ACDM2-B2 
(P = 0.1) models fit well the observational data at a relatively high 
significance level. 

We should conclude that the behaviour of the observed angular cor- 
relation function of the X-ray sources is sensitive to the different cos- 
mologies but at the same time there is a strong dependence on the bias 
models that we have considered in our analysis. By separating between 
low and high redshift regimes, we obtain results being consistent with 
the hierarchical clustering scenario, in which the AGN's are strongly 
biased at all cosmic epochs [10]. 



5. 



Conclusions 



We have studied the clustering properties of the X-ray sources us- 
ing the predicted angular correlation function for several cosmological 
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models. We parametrize the predictions for w(9) taking into account 
the behaviour of b(z) for a non-bias model, a galaxy conserving bias 
model with b(z) oc (1 + z) and for a galaxy merging bias model with 
b(z) oc (1 + z) 1 ' 8 . Utilising the measured angular correlation function, for 
faint and bright X-ray ROSAT sources, estimated in [20] and [1] respec- 
tively, we have compared them with the corresponding ones predicted 
in three cosmological models, namely the rCDM, ACDM and ACDM2 
(with a high cjg iass value). We find that the models that best reproduce 
the observational results are: (1) ACDM2 model {£l\ = 1 — £l = 0.7) 
with high o"g iass = 1.13 and bias evolution described by6(z)oc(l + z) 1 ' 8 
and (2) ACDM model (0 A = 1 - = 0.7) with low of ass = 0.9 and 
bias evolution by b(z) = 1. 
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